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O ' Abstract 

We study the sensitivity to the CP/T- violation search in the presence of ambi- 
guities of the theoretical parameters. Three generations of neutrinos are considered. 
The parameters whose ambiguities are considered are the differences of the squared 
■ masses, the mixing angles, and the density of matter. 

We first consider the statistics that are sensitive to the genuine CP-violation 
effect originating from the imaginary coupling. No ambiguity of the parameters is 
5_i ■ considered in this part. It is argued that the widely-adopted usual statistics are 



not necessarily sensitive to the genuine CP-violation effect. Two statistics that are 
sensitive to the imaginary coupling are proposed. The qualitative difference between 
these statistics and the usual one are discussed. 

Next we proceed to the case where the ambiguity of the parameters is present. 
The sensitivity of the CP-violation search is greatly spoiled when the baseline length 
is longer than about one thousand kilometers, which turns out to be due to the 
ambiguity of the matter effect. Thus the CP-violation search by use of CP conjugate 
channels turns out to require a low energy neutrino and short baseline length. It is 
also shown that such a loss of sensitivity is avoided by using T-conjugate oscillation 
channels. 
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1 Introduction 



The observation of the atmospheric neutrino anomaly by Super-Kamiokande ]]J provided 
us with convincing evidence that neutrinos have non-vanishing masses. There is another 
indication of neutrino masses and mixings by the solar neutrino deficit [0, ^, |, [5J, 

These results give us the allowed region and excluded region for the mixing angles and 
the mass square differences. Let us now assume that there are three flavors of neutrinos 
and denote the lepton mixing matrix U, which relates the flavor eigenstates u a (a = e, fi, r) 
with the mass eigenstates z/j with mass (i — 1, 2, 3) as u a — J2i=i U a i^i, by 



U 
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Here c^- and Sy stand for cos 6^- and sin 6^, respectively. The observations of the atmo- 
spheric neutrino anomaly give us an allowed region for sin #23 and the larger mass square 
difference (= Sm^). The solar neutrino deficit provides allowed regions for sin 6*12 and 
the smaller mass square difference (= Srri^). On the other hand, the no-oscillation results 
of reactor and accelerator experiments give us an exclusion region for sin #13 (e.g. ref. H). 
There is no constraint on the CP-violating phase 5. 

The idea of neutrino factories with muon storage rings was proposed to determine 
these mixing parameters (and the sign of Sm^ in addition). It attracted the interest 
of many physicists fTO, 0, p], O, HI |i~5|, HI HI HE HI |2( 
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and the neutrino factories 
turned out to be a very promising candidate for the next generation neutrino oscillation 
experiments. We will be able to observe neutrino oscillations even if sin^ is as small 
as 0.01. We will also be able to detect the CP-violation effects in such experiments 
The possibility to observe CP violation through long baseline neutrino oscillation 
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experiments was discussed in Refs. [p3j |24| , |25| , Pq , \F7[ [28[ and many papers followed these 
works. 

By what observation can we insist that we measure the CP violation? CP violation is 
characterized by the intrinsic imaginary part of a coupling in a Lagrangian. The presence 
of an imaginary part of a coupling gives different properties to particles and antiparticles, 
and it is observed to be CP/T- violation effects. Hence we need to discuss CP/T violation 
using a quantity which is sensitive to the imaginary part of a coupling. We have to be very 
careful to construct such a quantity since there is an indirect sensitivity to the CP violating 
phase (which constitute not an imaginary part of couplings but a real part) through a 
unitarity |29| . As for lepton CP/T violation in the long baseline neutrino oscillation 
experiments, one of such quantities is the difference of event rate between CP/T conjugate 



channels. We must take care of matter effect |3(| for CP conjugate channels since it gives 
difference to the event rate too. Therefore we must take into consideration ambiguities 
of the parameters as the matter effect can mimic the genuine CP violation partially. We 
will show that the ambiguities of parameters spoil the experimental sensitivity.^] 

: It is also important to consider experimental systematic errors and backgrounds, but we do not 
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We formulate the treatment of the ambiguity of parameters within the statistical 
method. It is important to build a statistical quantity that is sensitive to the imaginary 
part of coupling. We propose a proper statistics, and show explicitly that the sensitiv- 
ity to CP violation effect changes by taking the ambiguities into account. To this end 
we estimate how large exposure (proportional to number of muons and detector size) is 
required to observe CP- violation effect by the neutrino factory experiment. The optimal 
experimental setup (muon energy and baseline length L) is shown through such con- 
siderations assuming three generation of neutrinos which account for the solar neutrino 
deficit and the atmospheric neutrino anomaly. 

This paper is organized as follows. In section 2 we consider statistical quantities 
which are proper to search for a genuine CP- violation effect. There we assume that the 
parameters such as %'s and Sm^s are known without ambiguities. We will discuss in 
section 3 the case where the ambiguities of the parameters are taken into account. We 
present the requirement on the number of muons and the mass of a detector to observe 
genuine CP-violation effect through measurements of CP-conjugate oscillation channels. 
It is shown that the introduction of the ambiguities of parameters greatly change the 
sensitivity to the CP-violation search. In section 4 we investigate CP-violation search 
using T-conjugate oscillation channels. Ambiguities of parameters are taken into account 
also in this section. The sensitivity in this case is far better compared to the case using CP- 
conjugate channels, showing that this case is ideal (preferable) to search for CP- violation 
effect. Finally a summary and the discussion is given in section 5. 



2 Direct observation of CP-violation effect 

Let us first discuss the physical quantity which characterizes the presence of CP/T vi- 
olation in the long baseline neutrino oscillation experiments. Such quantities must be 
sensitive to the imaginary part of the coupling. We need to be particularly careful when 
we make use of CP-conjugate channels in the presence of the matter effect. 

Let us first recall how the imaginary part of the lepton coupling gets into the oscillation 
probabilities. We use the oscillation in the vacuum as a simplest example. The oscillation 
probability from u a to vp in the vacuum is given by 

P(u a -> Vf,;E,L) = J2 \Upi^ L,{2E) K^\ 

i 

= J^Re U^U^ cos ^^ + J^lmU:^ aj U pi U; jS m^^, (2) 

where E and L are the energy of neutrinos and the baseline length, respectively. The 

consider them in this paper. We assume that we can determine all the quantities such as particle energy. 
The only error taken into account is statistical ones. 
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unitarity of U leads to 



Im U* al U a2 Up X U; 2 = Im U* a2 U a3 Up 2 U*p 3 = Im U* al U aZ U pz U* pi = J, (3) 
which allows us to write 

P(u a -> up; E,L) = J2 ^ U^U^UpiU^ cos + J^T sin - (4) 

The Jarlskog parameter J defined by eq.(|3]) vanishes when all the L^-'s are real, and the 
second term of eq.(^) also vanishes. The existence of imaginary part of lepton coupling 
gives non-vanishing J, and thus we need to observe the quantity which is sensitive to J 
(including its sign) to search directly for CP- violation effect. Note that this statement is 
independent of the parametrization of U such as in eq.(|I|). 

Now we consider the quantity which is sensitive to the Jarlskog parameter J in presence 
of matter on the baseline. We put some assumptions for simplicity to discuss this point. 
Suppose that we have same initial energy spectra for neutrinos and antineutrinos. Also 
we suppose that the antineutrinos have about twice larger data size so that the expected 
event numbers for neutrinos and antineutrinos are equal in no-oscillation case. Then the 
oscillation event numbers of neutrinos N(u a — > vp) and that of antineutrinos N(u a — > i>p) 
are expected to be equal in vacuum if CP is conserved. The two event numbers are in 
practice different due to CP violation (if any) and matter on the baseline. The difference 
of the two event numbers, 

AN(6) = N(v a -> up; 6) - N(u a -> up; 6) (5) 

is intuitively sensitive to the genuine CP violation. Here the CP-violating angle 5 is 
explicitly written. This quantity does not vanish due to matter effect, even in absence of 
genuine CP violation. We thus consider 

AN{5)-AN(6 ), (6) 

where 6 = 5 G {0, it} corresponds to the CP conserving case. 
We stress here that the quantity 

N(u a -> up; 6) - N(u a -> up; 5 ) (7) 

is not necessarily sensitive to the genuine CP violation. To compare eqs.(|]) and (|^), let 
us consider the oscillation probability P[y a — > up;8), which is related to N(u a — > up; 6) 
roughly by 

E 3 

N(u a -> up; 5) ~ j^P^a -> v p \ 5). (8) 
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It can be shown in the high energy limit f|||31|| 
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J 
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sin 26*12 sin 26*23 sin 26*i3 cos 6 1 



13 j 



fl3 I 



sin 2 ^23 sin 2 2^i3, 



2v / 2G F n e E, 



(9) 



(10) 

(11) 
(12) 



and n e in eq.(|12|) is the average electron number density in the matter on the baseline. 
Recall that the Jarlskog parameter J defined by eq. (|3]) is expressed under the parametriza- 
tion eq.(ffl) by 



J = sin 2#i2 sin 26*23 sin 26*13 cos 6*13 sin 5, 
which is related to j in eq. QlCf) by 

J = j sin S. 



(13) 



(14) 



Thus the third term of eq.(0) is the contribution from genuine CP- violation effect. Note 
again that this statement is also independent of the parametrization. 
We obtain from eq. (0) 



P{vy. -> v e \ a) ~ P(y„ -> u e ; 6 , a) 



5m 2 X L 
AE 



5m 2 31 



j(cos5=F 1) 



I 1 ( aL 
~ 3 \AE 



(15) 



2 The limit is valid when E u > Srh^L/A, where Sr 
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y/ (Sm^ cos 2$i3 — a) 2 + (<5™|.i sin 2#i 3 ) 2 . 
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(— sign for <5 = and + for 5q = ir) and 



P(vn -> 5, a) - P[v^ -> z/ e ; 5, a) 
Pfz/,,, — > z/ n : -5, - 



u e ; 5, a) 



Sm^L 
AE 



(ah f 2 _ _ 5m|i 



^21 



j sin S 



+ o 



E 4 



-j cos 5 (2 cos 26*13 — 1) — 25 cos 26*13 sin 2 #12 



(16) 



We can observe in eq.fllS]) that the leading term of Piy^ — > v e \ 6, a) — Piy^ — > v e ] $o, 0) 
does not contain the genuine CP-violation term. This leading term can be canceled by 
taking the difference between the probabilities of neutrinos and antineutrinos. Equation 
(|16D is indeed sensitive to the j sin 5 term, though it contains an unavoidable matter effect 
term in addition. 

Our viewpoint is that the CP- violation search must be carried out directly by observing 
the contribution of the j sin 5 term in eq.(|D|), which originates from the imaginary part of 
the coupling as in eqs.(|3]), ([13]) and fll4|). This term is not the leading term at least in high 
energy region, but it can be picked up as a leading term by taking the difference between 
neutrinos and antineutrinos as in eq.([Tj]). Applying this consideration also to event rates, 
we regard that eq.@ is a better quantity than eq.(0) to pursue the possibility of direct 
CP- violation search. An analysis using eq.(|7]) is an usual parameter fitting method. It 
does not take into careful consideration whether or not N(u a — > vp\ 5) is sensitive to 
imaginary part of the coupling. Even if one obtains a result that N(5) 7^ N(5q) in an 
experiment which is sensitive only to the real part, there remains a possibility to build a 
certain Lagrangian with totally real coupling which can reproduce N(5). In this respect 
N(5) 7^ N(5q) cannot be the definite clue of the presence of CP violation. 

Let us further exemplify the difference between eqs.(^) and (|7|). We consider a following 
toy setup of an experiment. A source of neutrino beam is muons which decay into 
neutrinos at a muon ring. The neutrinos extracted from the ring are detected at a detector 
if their energy E u is larger than a threshold energy E^. The detector has mass Mdetector 
and contains iV target target atoms, which are related as 



AWt = 6-02 x 10 34 Mdct f° r . 

target [lOOkt] 



(17) 



We assume the neutrino-nucleon cross section a is proportional to neutrino energy as 



a = a E u , 
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where 



00 



0.67 x 10 38 cm 2 /GeV for neutrinos, 
0.34 x 10~ 38 cm 2 /GeV for antineutrinos. 



(19) 



The expected number of appearance events in the energy bin < E v < Ej (j 

1,2,... ,n) is then given by 



Nj(v a 



7rm 2 , L 2 



E v f Ua {E v )P{v a ^vp;S) 



d£„ 



(20) 



where m M is the muon mass, and f Va (E v ) is the neutrino flux that is concretely given by 
eqs.(|68|) and (K9|). We define 



C 



and 



Rj(u a 



0" ^target 

7rm 2 M detector 



E v f Va {E v )P{v a ^vp;E v ,8) 



(21) 



2-1 



(22) 



so that 



1.14 
0.58 



£ 2 

1/^5) = N^M detector —CRj(v Q 



vp\8) 



x 10' 



iV„. 



M detector / E^/lGeV] 
[10 21 ]' [lOOkt] VL/ [1000km] 



xRjfa^vfiS) (23) 



(1.14 for neutrinos and 0.58 for antineutrinos). A quantity iV^Mdetector is normalized in 
unit of [10 21 • lOOkt] in eq. (|23|) . The value of unity in this unit is a quite optimistic one 
compared to the presently discussed values The requirement on Mdetector in this 
normalization must be about unity or less so that we can observe the CP-violation effect 
experimentally. 

The widely adopted \ 2 statistical quantity based on eq.([7|) is defined by 



£ 

3=1 



[Nj(5) - Nj(8 )Y 



3=1 



N 3 (5) ^ Nj(S) 

3=1 



Mdetector ~JJ C <j 



Rj(S) 



+ £ 

3=1 



[Rj(5) - RjjSo)]' 



(24) 



where 



(25) 
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R 3 (6) = Rj(u a -+Uf,;5), Rj (6) = Rj {v a -> %; 5) (26) 
and n is the number of bins. Similarly we define based on eq. (P) asQ 

xM = ^ m + _ * w i^ c ^ ^w+w ■ (27) 

where 

ARM^R^-R^S). (28) 
We need both 5^0 and 5 ^ ir to ascertain that CP violation is present. We thus define 

X\ = min xl(6o), (29) 

<5oe{0,7r} 

Xl = min xl(S ). (30) 

<5 e{0,7r} 



and require 



Xl > X^ %(2n), (31) 
X2 > Xl 0% (n) (32) 



to claim that the CP-violation effect is observable at 90% confidence level in the method 
with n energy bins; more details on statistics is found in the Appendix. Equations ( |3~1~D 
and fl32"D is equivalent to 



-/V^Mdetector > {N ^ Mdetector )min;l,90% 

_ l^ 2 X90%( 2 ^) 



(33) 



mm 



and 



iV^Mdetector > ( -^Md e tcctor)min;2,90% 

= i^ 2 . X 2 90% (n) 

CE l min (y^ [A^W-A^(tf 



(34) 



respectively. 

We present example plots of (iV At M d etector)min;i,90% and (A^M de tector)min;2,90% in Figs.0 
and |||. We adopt only single-bin method here because to divide the energy region into 

3 This quantity depends on a certain model through subtraction of the matter effect. However, if we 
can observe this asymmetry significantly, then we would be able to conclude that there is a genuine 
CP-violation effect in the real Lagrangian, even if the real theory is not the model that we assume. 
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6 11 16 21 26 31 36 41 46 51 56 61 



^/[GeV] 

Figure 1: A contour plot of the required data size to observe the genuine CP- violation 
effect. A quantity (A^M detector ) min;li9 o% defined in eq.(??) is plotted in unit of [10 21 -100kt] 
as a function of muon energy and baseline length. Smaller value of this value means the 
higher sensitivity. Here E t ^ = lGeV, and the case of 5 = ir/2 is presented. Other param- 
eters are taken as shown in eqs.([35|) and (^). The smaller value of (iV At M detector ) mill . li 9o% 
means the higher sensitivity of the CP-violation search. The use of xl leads to the higher 
sensitivity as gets larger. 
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250 



6 11 16 21 26 31 36 41 46 51 56 61 

I [GeV] 



Figure 2: A contour plot of (A r Al Md e tector)min;2,9o% in unit of [10 21 • lOOkt]. The parameters 
are the same as in Fig.[l]. A larger data sample is needed compared to Fig.[l|. Optimum 
muon energy and baseline length makes a sweet spot in the graph. 

some bins does not make the sensitivity better, especially in the case without ambiguities 
of the parameters. We will explain the reason in detail at Section |3~2j The parameters are 
taken as follows so that they are consistent to the present experimental limit: 

sin #i3 = 0.1, sin ^23 = — =, sin# 12 = 0.5, (35) 



Srr4i = 3 x 10" d eV 2 , 5m 2 21 = 1 x 10" 4 eV 2 . (36) 

CP- violating angle 5 is taken to be tt/2. Matter effect a is approximated to be constant 
on the baseline, but its value depends on the baseline length since the longer baseline gets 
deeper in the earth. The Preliminary Reference Earth Model [3~2~| is adopted to estimate 
matter density as in Fig.|| |3l] . 



There is qualitative difference between Figs.|T] and|2|. Figure [l] shows that the sensitivity 
in terms of N(8) — N({0, tt}) enhances as the muon energy gets larger. This does not hold 
for AN (6) — AN({0,7r}) as seen in Fig.^. There is a sweet spot in this case in terms of 
muon energy and baseline length that optimizes the sensitivity to the CP- violation effect. 

We discuss the relation between xf and x\- In the CP/T- violation search, one compare 
(A^-(5),A>-(5)) with (iYj-(<Jo),#j(tfo)). One can equivalently compare (Nf sbal (6),AN j (6)) 
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Figure 3: Approximated matter density as a function of baseline length, calculated from 
the Preliminary Reference Earth Model |32| . 



with (iV Total (5 ), AJVj(<J )), where iV Total (5) = Nj{S) + Nj{6). A X 2 statistics defined by 



XA°o) = 2^ Nj otal {5) 



[AN 3 (6) - AN 3 (6 ) 
Nj ota \5) 



(37) 



obviously corresponds to Xi(^o)- The second term of eq.(^) is xl(^o) itself. Hence 
we focus to the first term to understand the relation between \\ an d X% We note that 
iV Total (5) is insensitive to CP-violation effect, since the magnitude of genuine CP-violation 
effect for neutrinos and antineutrinos are identical with opposite sign. The term we are 
discussing is thus insensitive to the imaginary coupling. Our aim was a direct CP/T- 
violation search or a search for an imaginary coupling, and thus we dropped this term to 
obtain x 2 . On the other hand, the dependence on the energy of jV Total and AN in high 
energy region is given by 



T Total 



AN 



(38) 
(39) 



which follows from eqs.(§), (|15D and ([161) . Hence the first term gets larger as energy gets 
larger, and dominates the right hand side of eq.fl3"P|) in the high energy region; total fit 
gets better in spite of a poor fit of the imaginary coupling. The high sensitivity obtained 
by use of \\ an d shown in Fig.|TJ was achieved in this way. 

The sweet spot seen in Fig.|2] can be intuitively understood. The CP-violation effect 
appears when the number of generations is more than three . On one hand the heaviest 



state decouples from the oscillation at the low energy region such that Sm^L, and 

on the other hand the lighter two generations are effectively degenerate in high energy 
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region such that > Sm^L. Thus the suitable energy region to observe CP-violation 
effect is roughly given by Sm^L < < Sm^L. The sweet spot exactly lies in this region 
reflecting that \\ * s indeed sensitive to the imaginary coupling. 

We mention here that the sensitivity for the case 5 = 7r/2 and for 5 = —tt/2 is not very 
different, contrary to the discussion by other authors [HJ. The previous works compared 
N(8) with N(5 = 0) alone, but we compared N(8) with both N(5 = 0) and N(5 = vr). 
One should keep in mind that CP symmetry is conserved not only in the 8 = case but 
in the 5 = it case; the imaginary coupling is absent in both two cases. The real coupling 
is different for these two cases, and thus we need to distinguish an experimental result 
with both of them. We took these points into account by the definition eqs.(p9|) and ( j3Cf) . 
We present in Figs.f| and [5] the sensitivity plot similar to Figs.|l] and |2|, but this case for 
5 = — 7r/2. We observe indeed no qualitative difference between Figs.|I] and |] and between 
Figs.|| and |j, respectively. 

We have seen so far that we can extract the imaginary coupling by constructing the x 2 
statistical quantity as in eq. (p0|) . The construction was motivated by taking the difference 
of event rates of neutrinos and antineutrinos. We can build another quantity along this 
idea as 

N{5) - ^\n(5), (40) 



12 



250 



6 11 16 21 26 31 36 41 46 51 56 61 

#H / [GeV] 



Figure 5: Same as FigJ^, but here 6 = —n/2. 
which vanishes when 5 = 8q. A \ 2 statistics for this quantity is given by 



i 2 



2( , _^ [N(6 )N(6) - N(6 )N(6)\ 

jr[ N(5 ) 2 N(5) + N(8 ) 2 N(5) ' 1 ' 

The sensitivity condition with 90% confidence level in n-bin method is 

Xl > xlo%(n), (42) 



or 



-^Mdetectoj. > (iV /J Md etector ) m i n; 3 j 9o% 

_ lE l Xl 0% {n) 



C L 2 

mm 



R{6 )R(6) - R(6 )R{6) 



2 



5o£{tW 1^ R(8 ) 2 R(8) + R(5 ) 2 R(5) 



(43) 



Figure || shows an example plot of (A^M detector ) min; 3 j90 % for the parameters given by 
eqs.(|3"o]) and fl3T)|). The graph is similar to Figs.^| and [5], which are obtained from \\. The 
quantity x\ is thus also sensitive to the imaginary coupling and suitable as a statistics for 
the direct CP-violation search. 
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^/[GeV] 

Figure 6: A contour plot of (A r Al Md e tector)min;3,9o% in un it of [10 21 • lOOkt]. The parameters 
are the same as in Fig.[l| A sweet spot is seen in this figure as was also seen in Figs.^j and 



3 CP-violation search in presence of ambiguities of 
the parameters 

3.1 Sensitivity to CP-violation effect in presence of ambiguities 
of the parameters 

We used N(5o) and N(So) in the definitions of \\ an d xl- Exact values of mixing angles 
and 5m 2 's are required to obtain iV(<5 ) and N(Sq), but they are not known in practice. 
We discuss how the ambiguities of parameters spoil the sensitivity to the CP violation. 

The ambiguity of parameters is especially important when we make use of CP con- 
jugate oscillation channels. It is because the genuine CP-violation effect in this case is 
contaminated by the matter effect. An estimation of matter effect is required, and the 
ambiguity of parameters are obstacles to the estimation. For the better understanding, 
let us get back to Piy^ — > v c ) — Pipp, — > v B ) in high energy region given by eq.flH]). It 
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consists of two parts, matter effect part APMatter and CP- violation effect part APcpv ; 



APMatter = 2 



8m\ 1 L\ 3 aL 



2 „ 5m%i 

-Bcos26 13 + -^ 



AE ) AE 
1 



—j cos 5 (2 cos 2^i3 — 1) — IB cos 2#i3 sin 2 12 



(44) 



5m 2 .^L\ 3 <5mn 



APePvS _ 2 ^j J-^, (45) 

There is an ambiguity in APMatter due to the ambiguities in 8m 2 's, 6's and a. A sensitivity 
to CP-violation part AP CPV is lost if the ambiguity of AP Ma tter is larger than AP CPV 
itself. Ambiguity of all the parameters contribute to the ambiguity of APMatter- It is thus 
important to take into account ambiguities of all the parameters. It is expected that the 
ambiguity of APMatter is large when APMatter itself is large. Recalling that APMatter is 
proportional to baseline length L (due to the factor aL/(AE)), one should obtain a poor 
sensitivity in the long baseline region. It is important to consider the sensitivity to the 
CP-violation effect when the parameters are not precisely known. 

We improve the discussion given in the previous section and formulate how to take 
ambiguities of parameters into account in estimating the sensitivity to the CP-violation 
effect. Suppose that one uses the parameters {xi} = {#12, #23, #13, Sfh 2 ^, Srh^, a}, which 
are different from the true values {xj} = {#12, #23, #13, $ m 2i, & m \ii a }> to calculate Nj(8 ) 
and Nj(5o). One will estimate 

NjiSo) = Nj(u a -> m {xi}, S ) (46) 

instead of Nj(So), where Nj(u a — > isp; {xi}, So) is evaluated from eq. (p0|) . Then matter 
effect can be overestimated and hence the sensitivity to the CP violation can be spoiled. 

First we see this using x\- I n this case one will estimate the fake CP violation due to 
matter as follows, 

ANj(S ) = Nj(u a -> v p \ {xi}, S ) - Nj(9 a -> v p \ {xj, S ), (47) 
instead of ANj(So). We then obtain 
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*«»> = jL m + m («) 

instead of xi(^o) defined in eq.(P7[). Adjusting the parameters {x^ within the ambiguities, 
one can provide a better fit to the expected values in no CP-violation case by minimizing 
X 2 ■ One can nevertheless infer CP violation is present in 90% confidence level if one 
cannot make the value of x 2 smaller than xlo%- We thus generalize eq.(^) and define 

f 2 = mm xl(5 ). (49) 

8 £{0,ir};{x t } 
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Figure 7: A contour plot of (iV Ai M detector ) min; 2 i amb,90% i n un it of [10 21 • lOOkt]. The parame- 
ters are the same as in FigfJ. The required data size is much larger than the no-ambiguity 
case shown in Fig.^. The sensitivity is rapidly lost when the baseline length gets longer 
than about one thousand kilometers. 

A criterion that CP-violation effect is observable for 90% confidence level in n-bin method 
is given similar to eq.(^) as 

xl > xlo%(n), (50) 
which can be rewritten in terms of N^M detector as 
A^M dctcctor > (A^M detector ) min;2jambi90% 

_ 1 L 2 xlo%{n) 



CE l m . n (ys [AR j ({x i },5)-AR j ({x i },5 )} 2 



(51) 



<5oe{o,7r};{^} I jr[ Rj({xi}, 5) + Rj({xi}, 5) 



We present (iV M Md e tector)min;2,amb,9o% m Fig.[7| to observe the CP-violation effect in 90% 
confidence level. All the parameters {x{\ are assumed to have ambiguities of 10%, and 
their central value are taken as in eqs.(|35|), and Fig.R, so that 



0.09 < sin #i3 < 0.11, 

0.9 . ~ 1.1 

7! < sm 023 < 7!' 

0.45 < sin < 0.55, 

2.7 x 10~ 3 eV 2 < Sm 2 31 < 3.3 x 10" 3 eV 2 , (52) 

0.9 x 10~ 4 eV 2 < 5m 2 1 < 1.1 x 10~ 4 eV 2 , 

0.9a < a < 1.1a. 
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It is seen in both figures that genuine CP-violation effect is difficult to be observed when 
the baseline length is longer than about one thousand kilometers. An estimation of matter 
effect is obscured by the ambiguity when the baseline length is long, and CP-violation 
effect cannot be separated from matter effect. This result can be understood qualitatively 
by the following rough estimation. CP-violation effect is hidden by the ambiguity of 
matter effect when matter effect is large enough. We require 



AR 



CPV 



Matter 



< 1 (53) 



as a rough estimation to observe CP- violation effect.^ Putting eqs. (f44l) and (f|^) into 
eq. fl53|) , one obtain a condition on L as 

L < 4£„ 3(6m 2 21 /6m 2 31 )jsm5 

~ a 2 sin 2 6 23 sin 2 26 13 cos 26 13 ' 1 ' 

Applying our test parameters eqs. (|35"D , (|36D and Fig|| to eq.(|54]), one obtains 



L < 1250km, (55) 

which is consistent to Fig.0. 

Next we illustrate using x 3 that the sensitivity to the CP violation is lost in presence 
of ambiguities of parameters. The correspondent of eq.(|51|) in this case is given by 



Mdetector > ( -^ju ^detector ) min;3,amb,90% 
1 X-90% 
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CL2 [Rj({^}^o)R,({xr},5) - Rjdx^^o)^^}^) 



mm 



-(56) 



5oe{0,7r};{£i} | j-f Rj ({x;}, 5 ) 2 Rj ({xt}, 5) +Rj({Xi}, 5 ) 2 Rj({Xi},5) 

We present a plot of (A^M detector ) min;3iamb>90% in Fig.|. The ambiguity of parameters 
makes the sensitivity worse also in this case, as we see by comparing Figs.^ and [|. It is 
seen, however, that the sensitivity shown in Fig|| is better than that shown in Fig.[7], which 
means that x 3 avoids the ambiguity of matter effect better than \\. One can understand 
the better sensitivity of x\ as a cancellation of ambiguity of sin # 13 when the high energy 
limit applies. The dominant parts of i?'s are given in the high limit by 

Rjdxi}^) = B(S + T) + V 

Rj({xi},8) = B(S-T)-V 

Rj({xi},5) = B(S + T) (57) 

Rjdx&S) = B(S-T), 



4 One should actually use as a denominator the ambiguity of APMattcr, not APMattcr itself. We 
tentatively use eq.(p3|) to give a rough estimation, however. 
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6 11 16 21 26 31 36 41 46 51 56 61 

^/[GeV] 

Figure 8: A contour plot of (A^M detec tor)min;3,amb,90% in unit of [10 21 • lOOkt]. The param- 
eters are taken to be the same as Fig.[| Note that the sensitivity is enhanced compared 
to Fig.0. 



where 

--<^(^) s (a)i 



Only the ambiguity of sin #13 is taken into account here, and thus B = sin #23 sin 2 2^13. 
Using eqs.(|5TD, we obtain 

, (flj S in<S) 2 _ (jsin(5) 2 

Note here that B vanishes in eq . (|60|) . The ambiguity of sin# 13 is thus canceled away in 
the high energy limit. We also expect that the ambiguity does not spoil the sensitivity to 
the CP-violation effect even in the lower energy. 

The sensitivity can be enhanced by a construction of a good statistics such as xh but 
in general the sensitivity to the imaginary part of coupling is smaller as the baseline length 
becomes longer. We confirm it for a couple of parameter sets by presenting Figs.^j and [10 



It is also seen in Fig.|9| that the sensitivity for 5 = ir/6 and for 5 = 5ir/G is quite different. 
The genuine CP violation has a same magnitude for both of these two cases. On the other 
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hand, the term proportional to cos 5, which is contained in the matter effect term (see 
e.g. eq.(|44])), has an opposite sign. The magnitude of matter effect contamination is thus 
different, and it leads to the difference of the sensitivity according to our discussion that 
the sensitivity to CP- violation effect is controlled by the magnitude of matter effect. 

We have been keeping on discussing the direct observation of CP- violating effect. One 
can verify our results by use of a quantity similar to x\ defined in eq. fl24|) . Equation ( [24]) 
focuses on the difference between N(v a — > Vp;8) and N(u a — > isp;5 ), and xl nas little 
sensitivity to the genuine CP- violation effect as a result. Instead we define 

Xlasym(^o) 



2 



- N,M detectoI -C ^ m + ^ m j, (61) 

which is sensitive to the genuine CP- violation effect |34|].(] A quantity analogous to 
(A^M detector ) min;liamb9 o% is defined by 

(iV /i il^detector)min;lasym,amb,90% 

= X^ 0% (2n) 

CE l . {^[R j ({x i },5)-R j ({x i },-8)\* " [^({xi}^)-^-^},-^] 2 



mm 



{£} | ^ ^({xi},^) ^ ^ ( {x, } , 5) 



(62) 



Figure |TT] shows a contour plot of (A^M de t ec tor)mm;iasym,amb,90%- It is seen that the sen- 
sitivity is lost when the baseline length is longer than about two thousand kilometers, 
which is qualitatively consistent to the results obtained in this section. 



3.2 Energy dependence 

We discuss the binning of the neutrino energy in a search of CP violation. We recall that 
the genuine CP violation in terms of oscillation probability is given by 

5m%L 5m\ L 8m%»L , , 

J sin - 21 sin sin - i? . 63 

AE U AE V AE U y 1 

Applying Sm^i <C dm^ and 5m\ x L/ (AE V ) <C 1, eq.flBB]) is rewritten to be 



-J^sin 2 ^. (64) 
AE V AE V 1 ; 



5 This quantity requires both CP- and T-conjugate channels; we consider this quantity just to verify 
the discussions so far. 
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Figure 9: Same as Fig.|8], but for different parameters. All the graphs presented here are 
for (JmJi = 1 x 10 _4 eV 2 . The graphs in left column are for sin# 13 = 0.1 while the ones in 
right column are for sin# 13 = 0.05. The top two graphs are for S = 7r/6, the second two 
graphs are for S = 7r/2, and the bottom two graphs are for 5 = 5it/2. Parameters not 
presented here are taken to be same as Fig.|2|. The difference of the sensitivity for 5 = tt/6 
and for 5 = 5n/6 is due to the difference of matter effect. 
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Figure 10: Same as Fig.O, but for brr? 21 = 5 x 10 5 eV 2 . 
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Figure 11: A contour plot for (A^M detector ) min; i asymiambi9 o% defined by eq.(|62|). The sensi- 
tivity is lost also when the baseline length is longer than about 2,750 km. 

It can be seen from eq.fl64|) that the genuine CP- violation effect has a definite sign as a 
function of E v . It is also applicable to the event rate N(5). Dividing the event rates into 
energy bins is thus meaningless and unnecessary to search for the CP- violation effect, 
when the matter effect is absent. All one need to do is to observe the total counts of 
neutrinos. This is of practical importance for experimental studies since the determination 
of neutrino energy is very challenging . 

On the other hand, a single-bin analysis does not necessarily remain advantageous once 
the matter effect is taken into account; the multi-bin analysis is required to remove the 
matter in such a case that the considered \ 2 is sensitive not only to genuine CP-violation 
effect but to matter effect (see Fig.[T2|) . The number of events per bin is sacrificed by 
the bin dividing, and each bin has a relatively small number of events compared to the 
single-bin analysis. As a result, the best-fit point of multi-bin analysis is less robust than 
the single-bin analysis, i.e., the best fit point of multi-bin analysis easily lies far away 
from the true parameter point. 

We conclude as follows from the above considerations. Experiments to search for CP- 
violation effect should be made with the setup where the single-bin analysis gives the best 
sensitivity, which means that the matter effect contamination is small. Such a setup has 
another practical advantage in addition: measurements of the neutrino energy are not 
required in single-bin analyses. One need not take care of the correlation between bins 
due to finite energy resolution, which makes the sensitivity to the CP violation in the 
experiment worse. 
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Figure 12: The most effective binning method to observe the CP- violation effect with X3. 
This corresponds to Fig.[8[ Single-, 3-, 5-bin analyses are compared. White, light gray, and 
dark gray regions mean where single-, 3-, 5-bin analyses are most effective, respectively. 
This shows that multi-bin analysis is required in long baseline range. 



3.3 Dependence on sin #13 of the sensitivity to the CP-violation 
effect 

We finally discuss the correlation between parameters on the sensitivity to the direct 
CP-violation search. 

The magnitude of CP violation is determined a single parameter, namely the Jarlskog 
parameter jsin5. The correlation between S and other parameters such as #13 is often 
discussed, but it is heavily dependent on the parametrization. The presence or absence 
of CP/T violation can be determined without any correlations to the mixing angles.^ 

We present Figs.O and in] to show the correlation between sin # 13 and 5. They are 
sensitivity plots using where 5 and sin #13 are varied while E and L are fixed. Figure 
is a test plot for E = lOGeV and L = 1000km, and Fig.|TJ is for E = 20GeV and 
L = 2000km. A direct CP-violation search is expected to be possible with this setup, 
as seen in Figs|| and 10. It is expected in these figures that the sensitivity scarcely 
depends on # 13 if the statistics is correctly sensitive to the genuine CP- violation effect. It 



6 The real part of coupling is in fact another intrinsic parameter which is independent of the parame- 
terization. The correlation between real part of the coupling and imaginary part of that will be present. 
This is the only possible correlation for the CP-violation effect. 
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Figure 13: A similar plot to Fig.[7], but this time = lOGeV and L = 1000km are fixed 
while sin # 13 and 5 are varied. The contour is nearly vertical, which reflects the fact that 
the value of sin # 13 is not important to consider the sensitivity to CP- violation searches. 



is illustrated by a rough estimation of xl m such 



A 3 



(AP< 



CPVJ 



L 2 [{8m 2 x / 5m 2 x )j sin 8\ 



E 
L 2 



B 



5m 2 21 



E \ Sm^ 



[sin 26*12 sin 2#2 3 sin 29 13 cos 9i 3 ) 2 



sin #23 sin 2 29 



13 



E V Sm^ 



,)m21 ] sin 2 29 12 cos 2 9 23 cos 2 9 v> . 



(65) 



Equation ( |55D depends on 13 only through cos 6*13, which is almost unity for 9 13 <^ 1. 
The dependence on 13 should be thus quite small if the statistics is sensitive to jsinS. 
We compare Fig.|13| with Fig.|TJ| to confirm the above discussion. It can be seen in Fig.|8| 
that the parameter for Fig.[13| is more sensitive than that for Fig.[T3. We see that the 



dependence of the sensitivity upon sin #13 in Fig. 13 is quite small, while the dependence is 



larger in Fig.[T4|. This is indeed consistent to the above discussion; the larger matter effect 
gives a sizable contribution to the numerator of \\ in the latter case, and the estimation 
given in eq.(|65|) does not apply We thus conclude that L ~ 1000km and E ~ lOGeV is 
the optimum setup to search for a direct CP-violation search with use of the statistics 
given by eq. (|4i"|) .f\ 

7 An optimum setup should change if one can find other better statistics, since the sensitivity itself 
depends on the adopted statistics. The difference of Figs.0 and || is an example. 
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Figure 14: A similar plot to FigO, but this time E„ = 20GeV and L = 2000km. Matter 



effect is larger in this parameter compared to Fig.|H| which leads to the dependence of 
the sensitivity upon sin# 13 . 



4 T- violation search in presence of ambiguities of pa- 
rameters 

We have discussed in the previous section that the ambiguity of matter effect spoils the 
sensitivity to the CP- violation effect. One can then expect that one can avoid the loss of 
sensitivity by use of T-conjugate oscillation channels, which is free from the matter effect 
J55 . 

It is convenient to redefine 



and 



Nj({xi}, 5) = N{v a -> v p \ {xi}, 5) 



Nj({xi}, 5) = N(vp -> u a ] {x^, 5), 



(66) 



(67) 



where a,/3 = e,/i, r, and v denotes neutrinos and antineutrinos collectively. We are to 
consider initial neutrinos of different flavors, and we must take into account the difference 
of the energy spectra of the neutrino beam. For example, the flux of u e and obtained 
from the decay of muons with energy is given in terms of x = E v j 'E^ by 



f Uc (x) = 12x 2 (l-x) 



(68) 



and 



f u (x)=2x 2 (3-2x), 



(69) 
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Figure 15: A contour plot of the required data size to observe the T- violation effect. A 
quantity (A^M detector ) min;Tiambi9 o% defined in eq.flTDp is plotted in unit of [10 21 • lOOkt]. 
The parameters are the same as in Fig.^|. This figure is similar to Fig|| and it is seen 
that the use of T-conjugate channels is robust to the ambiguities of parameters. 



respectively. The quantity xi defined by eq.([|l|) is suitable in such a case. We define 
(A^Mdetector)min;T,amb,9o% i n the same manner as the rightmost side of eq.flSB]) so that 



detector Jmin;T,amb,90% 



1 "^M ^90% 



n 



CL2 ^ /A [i2 i ({5 i },5o)^({^} ) 5)- J R i ({^},5 )i2 j ({^},5) 



,(70) 



mm 

5 e{0,7r};{iJ Rj{{Xi},^) 2 Rj{{Xi},S) + J R i ({xJ,(5o) 2 J Rj({x i },(5) 



while this time Rj(5) = Rj{y a — > 8) and Rj(S) = Rjivp — > v a \ 5). T- violation effect is 
considered to be observable when 

A^Mdetector > (^^detector) min;Tiamb 90% (71) 

is satisfied. 

We present in Fig.|l5] a test plot of (A r A1 Md e tector)min ; T,amb,9o%- The parameters are 
taken as shown in eq.(|5*2]). It can be seen that Fig.|T3| is qualitatively similar to a plot in 



absence of the ambiguity of parameters, which is presented in Fig.|[ A sweet spot, which 
is expected from the naive estimation in terms of oscillation probability, still remains in 
Fig.|15|; thus we find that the CP/T- violation search via T-conjugate channels is robust 
to the ambiguity of the parameters. 

Longer baseline length is in general preferable when T-conjugate channels are available, 
since eq. (|65[) applies without being troubled with matter effect contamination. This is 
in contrast to the CP-conjugate case, where the matter effect obscures the genuine CP- 
violation effect when the baseline gets too long. 
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5 Summary and Discussion 



We discussed the optimum experimental setup and the optimum analysis to see CP- 
violation effect. 

We first discussed the difference between the direct measurement and the indirect 



measurement [ 29] . Genuine CP- violation effect is characterized by an imaginary part of 
couplings in Lagrangian and hence quantities sensitive to this imaginary part should be 
used to measure the CP violation. To see this we introduced two statistical quantities, 
Xi (eq.@) and X \ (eq.(U)). 

Usually xi is used in analyses of neutrino factories. We can test using this whether 
the data can be explained by the hypothetical data calculated assuming no CP-violation 
effect. We saw, however, that this quantity is not necessarily sensitive to the CP- violation 
part of the coupling. In high energy region it is sensitive almost only to the CP conserved 
part of the oscillation probability. We can tell about genuine CP-violation effect only 
thorough unitarity relation of three generation. The sensitivity to the CP- violation effect 
is often indirect. Thus we concluded that we should not use it to measure CP violation 
since it often measures the CP-violation effect only indirectly. 

On the other hand, we can test with x\ whether the asymmetry of oscillation prob- 
abilities of neutrinos and antineutrinos exists. We have seen that x\ is sensitive to the 
CP violating part of the oscillation probability, and thus it is more suitable quantity to 
measure the CP- violation effect. 

We saw the relation between x\ and \\ and found that to pick up an imaginary part 
of couplings we need to see the difference between particle and antiparticle. In this sense 
we also introduced a statistics x\ (eq.(f4TD). This statistics gives better sensitivity to 
measure CP-violation effect directly when we consider the ambiguities of the theoretical 
parameters. 

Then we investigated the influence of the ambiguities of the theoretical parameters 
on \\ and x\- Since the matter effect causes between the difference the oscillation prob- 
abilities neutrinos and antineutrinos, we have to estimate fake asymmetry to search for 
the CP- violation effect. However, we will always "overestimate" the fake CP violation 
because of the ambiguity of the theoretical parameters, and hence we will always esti- 
mate the genuine CP- violation effect too small. The matter effect increases as baseline 
length increases, and we will lose the sensitivity to the asymmetry due to the genuine 
CP-violation effect in longer baseline such as several thousand km. 

The sensitivity of x\ to genuine CP violation is lost much more than that of xl- This is 
due to the partial cancellation of the ambiguity by sin #13. The ambiguity of estimation of 
matter effect is partially canceled in the numerator. We found that x\ is better statistics 
to see CP-violation effect directly. 

Using xi we studied the correlation between sin 13 and 5. Comparing Figures [13] and 



IJ, we found that in general that we have better sensitivity to CP violation with baseline 
length 1000km than 2000km. Moreover, if the statistics is only sensitive to the imaginary 
part of couplings, Jarlskog parameter J, it has no dependence on sin 6*13.0 Indeed in Fig.|T3 



3 As long as sin $13 term has dominant contribution. 
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we see this behavior while in Fig.|l4| we see strong dependence on sin #13 of the sensitivity. 
In this sense we also understand that baseline length 1000km is better to see CP violation 
directly. Furthermore, Fig|| and Fig. 10 show that around 1000km is optimal baseline 
length for various parameter sets. 

Taking the statistics which is sensitive to the imaginary part of the lepton couplings, 
we first showed that there is a sweet spot in terms of and L when the ambiguities of the 
parameters are not considered. We have then taken the ambiguities of all the parameters 
to be 10% and showed that the sweet spot survives in such a case. We expect that the 
sweet spot also survives when we adopt the more realistic values of the ambiguities. We 
optimistically expect that other parameters will be determined with ambiguities less than 
10% except for #13 in the future. The large ambiguity of #13 is seemingly enough to wash 
out the sweet spot. We have mentioned in Section |3.1j , however, that the ambiguity of #13 
is canceled by use of the statistics x|. We thus conclude that the experimental setup of 

~ lOGeV and L ~ 1000km is desirable even in the real experiment. 

We finally studied T asymmetry. There is no fake asymmetry due to environmental 
effect such as the matter effect. We found that the naive expectation on CP- violation 
phenomena is indeed realized. 

It is required to find another way to see CP-violation effect if we can observe only 
appearance events of v e — > and v e — > v^. Otherwise we cannot observe CP- violation 
effect in neutrino factories with long baseline (> 1000km) as the asymmetry between 
neutrinos and antineutrinos. On the other hand, we can observe CP-violation effect as 
the T asymmetry very well. Therefore it is very important to establish a way to observe 
this asymmetry experimentally. 



Appendix: Statistics 

We explain a detail of the statistics used in this paper to estimate how many events we 
need to tell the existence of the genuine CP- violation effect. To state the feasibility of the 
experiment we consider not only how well we can distinguish two theories (two parameter 
sets) but also how well the best fit point lie in the true value. For example, even if in 
nature 5 = tt/2 is realized, we are not sure that the best fit point for S sit there. We will 
have the best fit point value other than 5 = tt/2 and hence we have to take care of this 
possibility to state the feasibility of the experiment. 

To estimate it, we employed the concept of the "power of test". In the test, we set up 
"null hypothesis", H , which should be rejected and its "alternative hypothesis" against 
H , Hi. In this paper we were interested in whether we can insist the existence of the 
CP-violation effect, and hence we set the null hypothesis, 

H : 6 = 5 (72) 

and "alternative hypothesis" against H , 

Ht : 8^S . (73) 
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We also define "test statistics" to give a criterion to reject H for a real data set Nf x . 
In this paper we examined three test statistics corresponding to \\iX\i an d x! : 

Tl{6o) = \ WW) + \ WW) ' 1 = *oX} Tl(5o) ' (74) 

T (5) = y [W! - g> - WW - ^ h (^°)}] 2 T — min T (5 ) (75) 



[iV? h (5o) x iVf x - N^(S ) x iVf x12 



3{o) = ^JWW)rWWWWm ^, 3 = ^X} T3(5o) ' (76) 

where N^ h is the event number assumed by the theory with S. Hereafter we use as 

an example \\ f° r the explanation. Furthermore for simplicity we abbreviate T\ (and 
accordingly \\) as 

T '-\ WW) ■ () 

To reject H at a "level of significance," we require 

Ti > Xl(2n). (78) 



Then the question is how well the inequality ([78] ) is satisfied for given value 5. This 
probability is called "power" ; 

P l (S) = P 5 (T 1 >j^(2n)). (79) 

This is the probability that we succeed in seeing the CP- violation effect in the experiment. 
Thus we have to require that this probability should be larger than 7, which is almost 1. 

To estimate the probability, often we generate event sets with a given event rate and 
check whether H is indeed rejected according to the inequality ( [78] ) with the probability 
7-0 Instead to do so, here we make a following approximation. First, we approximate Ti 
as 

Ti ~ \ ww) \ ww • (80) 

where iV* 1 ^^) is "the maximum likelihood estimator", i.e., 

» [ivr - jg^gg < A wr - ^ h (^)] 2 (RU 

V m h (S B ) "V ^(5) ' [ * } 



9 We have to generate enough event sets enough to conclude that Hq is rejected with the probability 

7- 
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for all 5. Equation © holds well if |iVf x - N^(6 B )\ < 0(^N} h (5 B )), i.e., the fit for the 
data iVf* by A^ th (5 B ) is good enough, and N^ h does not vary so rapidly around 5b- We 
also assume that the estimator is almost the true value, i.e., 5b — 5. Thus, 

- [ N r - Afffl] 2 a m\5) - jg (^o)] 2 

1 V iVf(5) Af(5) 

"2. ^ (82) 

With this approximation we calculate the power ( |T9"D as follows: 

aw = « (t |JVf, ^S (i)1 ' > *«< 2 "> - *i) ■ < 83 > 

The left hand side in parenthesis of P$ follows the x 2 distribution with 2n degrees of 
freedom so the requirement that the power f3\ (5) should be larger than 7 is equivalent to 
the condition 

X \>xl(?n)-xl l (?n-f) 1 (84) 

where / means the number of parameters included in the theory. For example, if we take 
the 0.1 level of the significance and the require the power to be 0.99 level, then 

Xi>x2.i(2n)-xg.99(2n-/)- (85) 

Since in general if 7 ~ 1 then Xj(^ n ~ f) is very small for small n, it is omitted in this 
paper.Q Thus we required^ 

Xl > Xl(2n). (86) 
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